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I. Introduction 

Aharonov and Bohm, in their historic 1959 paper jjj, 
showed theoretically that charged matter waves are phase 
shifted if they pass through a magnetic vector potential 
even when there is no magnetic field present in that re- 
gion of space. This so-called Aharonov-Bohm effect (AB 
effect) is a purely quantum mechanical phenomenon. The 
magnetic vector potential has no effect on classical parti- 
cles. The AB effect has been observed in several experi- 
ments using electron interferometryQ. 

The AB effect can also be studied in confined geome- 
tries with leads by looking at behaviors of electron trans- 
port through these systems. An one-dimensional struc- 
ture such as a small metal ring with a confined magnetic 
flux has been used to investigate the connection between 
the conductance oscillations and the AB effect Re- 
cent developments in micrometer-scale technology have 
also made possible new ways to observe the electron trans- 
port through confined two-dimensional cavities. These 
ballistic electron waveguides are formed from the two- 
dimensional electron gas (2 DEG) found at the interface 
of GaAs/AlGaAs heterostructure[Q. They are formed at 
very low-temperature and with sizes of less than /zm to 
make electrons phase-coherent inside. Since the quantum 
mechanics can be used to explain the behavior of elec- 
trons, these devices also provide a natural place to look 
for the AB effect 

In these systems, the magnetoresistance has been ob- 
served to oscillate with a period of h/e or h/2e as the 
magnetic flux is varied. In the ballistic regime, h/e os- 
cillation is expected due to gauge invariance in the quan- 
tum mechanics. In the diffusive regime (when both the 
elastic mean free path and the phase-coherence length 
are not larger than the size of the sample), oscillations 



with period h/2e were observed, and had been predicted 
theoretically [Q. Kawabata and Nakamurajl3[ argued 
that h/2e oscillation is theoretically predicted even in bal- 
listic cases using the semiclassical scattering theory. 

In this paper, we study the conductance properties of 
a two-dimensional ballistic circular billiard with two leads 
attached. Circular billiards have been studied before in 
several contexts. LinJTJ] calculated the conductance for 
this system using the semiclassical conductance formula. 
In Refs.(ll| and fl6|| , circular billiards with two leads were 
investigated with a homogeneous magnetic field applied. 
We will apply a point magnetic flux at the center of the 
circular billiard (AB billiard). Then the cavity contains 
only the vector potential with no magnetic field. In the 
system we consider here, the size of the cavity and the 
temperature are small enough to make the system ballis- 
tic, the boundary is formed by a hard-wall potential, and 
electron-electron interaction is ignored. Under these con- 
ditions, the Landauer formula ]l7|]-[|T9| f or two leads can 
be used to obtain the conductance of the system. (Be- 
cause h/e oscillation is expected in this ballistic model, 
we will only look at the magnetic flux between — h/2e and 
h/2e.) In addition to considering the effect of the vector 
potential on conductance, we will also consider the effect 
of the shape of the leads on conductance, and investigate 
conductance for two different shapes of leads, straight and 
conic (see Fig. 1). We will find that the shape of the leads 
and the relative position of leads can have an important 
effect on the conductance. We will also find that the AB 
effect can dominate the system under a special circum- 
stance. 

In Sec. II, we first consider the properties of eigen- 
states of a closed circular billiard with a magnetic flux 
line down its center. In Sec. Ill, the open billiard is stud- 
ied, and the boundary element method which is used to 
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obtain the conductance is discussed. In Sec. IV, we show 
our numerical results. Finally, we summarize in Sec. V. 

II. Closed Billiard 

In this section, we obtain the energy levels for an electron 
in a closed circular billiard with a magnetic flux line down 
the center. Reimann et al. |^0| calculated the energy den- 
sity of states semiclassically using the Gutzwiller's trace 
formula. In this paper, we must understand this system 
quantum mechanically in order to properly interpret the 
behavior of the open billiard with leads attached in the 
low energy range. In Ref. we studied this closed 

system when the flux line oscillates periodically in time. 
Then magnetic and electric field are induced in the cavity 
and the system can undergo a transition to chaos. When 
the flux line is constant in time, only a magnetic vector 
potential appears in the cavity and there is no chaos but 
we can study the AB effect. 

The Hamiltonian of the system can be written 
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where a = $/$o is a dimensionless quantity when $ is 
the magnitude of the magnetic flux through the center of 
the cavity and $o = h/e, and Vii(r) is a hard-wall po- 
tential confining the electron to the interior of a circular 
billiard of radius R. From here on, we will only consider 
a in the range —0.5 < a < 0.5, which is sufficient because 
of gauge invariance. For a — and a = 0.5, the energy 
eigenvalues are two-fold degenerate. However, the flux 
breaks the time-reversal symmetry of the system when 
< |a| < 0.5, and removes these two-fold degeneracies. 
After solving the time-independent Schrodinger equation, 
Ho(a)\l, n; a) = Ei n (a)\l, n; a), we can find energy eigen- 
states, 



(r,9\l,n;a) = J|;_ c 
and energy eigenvalues, 
Ein(a) 
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where m is the mass of the electron, and (Ji n (a) is the 
nth zero of J\i- a \ (x) (I: integer, n: positive integer). The 
energy eigenvalues, Ei n {a), are plotted in Fig. 2 as a func- 
tion of a. We can easily prove that the energies are sym- 
metric in a. We also see no level-repulsion as expected 
since this system has no chaos. When a = 0, |Z, n; 0) and 
| — I, n; 0) are degenerate, and when a — 0.5, \l, n\ 0.5) and 
\l — 1, n; 0.5) are again degenerate. 

III. Open Billiard 

We now create an open system by attaching two infinite 
leads to the circular AB billiard. As shown in Fig. 1, we 



will study two different lead geometries. For one geom- 
etry, the leads are straight. For the other geometry, the 
leads have conic shape. For both geometries, there are 
four dimensionless quantities which characterize the sys- 
tem. There are two geometric parameters, A (the opening 
angle of leads) and 7 (the relative angle between leads). 
There are two physical quantities, a = < E > /'3 ) o (the di- 
mensionless magnetic flux) and e — (mR 2 /h 2 )Ep (the 
dimensionless Fermi energy). 

For straight leads, we can calculate the number of 
propagating modes if we are given values for A and e. 
If we define the Fermi wave number as kp = \/2mEp/h 1 
there are n propagating modes if 
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where W — 2R sin(A/2). In other words, if e n < e < 
e n+ i, then there are n propagating modes in leads where 
e„ is defined as 

(5) 



8 sin 2 (| A)' 

All the remaining modes are evanescent modes. Here we 
are interested in the energy range, eo < e < £2- 

For conic leads, there is no way to distinguish be- 
tween propagating and evanescent modes before solving 
the problem, because all existing modes are propagating 
modes. But, after solving the equation, we find that, ex- 
cept for a finite number of lowest modes, all higher modes 
will be reflected completely. As we will see in our sub- 
sequent results, transmission is not necessarily zero when 
e < ei for conical leads, and the second propagating mode 
starts to transmit when e is lower than £2- We will call this 
"tunneling." The two lead geometries give similar results 
except for this tunneling effect. 

We will briefly look at the numerical method. The 
partial differential equation we are solving is a Helmholtz 
equation, 

(v 2 + fc!o*M) = 0, (6) 

with the boundary condition, ^lon ail boundaries = 0. The 
boundary element method |^2|, is used to solve this 
problem. We define a loop C consisting of two circu- 
lar segments, Pi and Pi, and two segments, C\ and C2 
across the contacts between the cavity and leads (see Fig. 
I). For straight leads, C\ and C2 are linear segments, 
but for conic leads they are parts of a circle, making C a 
complete circle. The loop C divides the whole region into 
three parts consisting of the interior of the circular cavity, 
the lead I, and the lead II. We will denote basis states for 
the interior of the circular cavity by {</>| c ^ : integer}; for 

lead I, {4>ni \ n i > 1}; an d for lead II, {0„ 2 |n 2 > I}; 
where + (— ) stands for an outgoing (incoming) mode. All 
basis states satisfy the Helmholtz equation. For example, 
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and so on. For both lead geometries, (j)^ is given by 

4>\ c \r,6)=c\ c) Jn- a \(k F r) e m . (8) 
For the straight leads, we have 
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where (a;j,yj)'s are local coordinates, k ni = 
^{2m/h 2 )E F - ( ni 2 TT 2 /W) 2 , and 



K ni = ^(ni 2 TT 2 /W) 2 - (2m/h 2 )E F . (In this case, a 
is set to zero inside the leads as an approximation to 
simplify the bases, i = 1,2.) We have N p propagat- 
ing modes where N p is the largest integer n satisfying 
(n 2 ir 2 /W) 2 < (2m/h 2 )E F . 



For the conic leads, we have 
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and 



(11) 

where (rj,^)'s are local coordinates and i?^ and iJ^ 2 ' 1 
are the Hankel functions (i = 1,2). The normalization 
constants, Cnj's, are chosen to make the incoming flux of 
electrons unity. 

The solution to the Helmholtz equation, Eq. (6), can 
be expressed as follows. Inside the cavity, 

oo 

*( C >M)= E «^ (C) M); (12) 

I— — oo 

in lead I with an incoming wave of kth mode, 
^(r 1 ,e 1 )=4 ) ~(r 1 ,6 1 ) + f2r jk ^ 1) + (r 1 ,e 1 )-, (13) 

and in lead II, 

oo 

* (2) (r 2 ,02) = E*^f (r 2 ,e 2 ), (14) 

where and are reflection and transmission ampli- 
tudes, respectively. 

Let us now consider the boundary conditions on loop 
C. (The boundary conditions on the walls of leads are 
already taken care of, because the bases we chose in leads 
satisfy those conditions.) On boundaries C\ and C*2, we 
have the boundary conditions, 
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where n is a normal coordinate of boundaries (i = 1,2). 
To express the boundary condition on Pi and Pi , we will 
introduce bases, : integer} and {£^1^2 : integer}. 

They are defined on Pi as 



^(A.) 



(17) 



where is a longitudinal coordinate along the boundary 
Pi and Aj is the size of the boundary Pi (i = 1, 2). Then 
these boundary conditions are written as 
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We multiply Eq. (7) by *< c ) and we multiply Eq. (6) 

(c) 

by cpl , and subtract one from another. Then, we inte- 
grate over the loop C. From the Green's theorem, we get 
the equation, 
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*< c ) and are obtained using Eqs. (15), (16), (18), 

and (19). With suitable truncations, after integrating 
over each segment, we finally obtain a matrix equation 
A • x = B. Solving this matrix equation gives us the val- 
ues for Tij , Uj , b$ , and b$ . (a; 's can be obtained using 
these values following similar steps.) 

The reflection and transmission amplitudes, r%j (i, j < 
N p ) and tij < N p ), respectively, determine the left 
half of the S matrix, 



t' 



(21) 



for the open circular billiard. In Eq. (21), r, t, r', and t' 
are (N p x A r p )-submatrices of S, and N p is the number of 
channels in the lead. The quantities r and t are matrices 
containing the reflection amplitudes and the transmission 
amplitudes, respectively, when there is an incoming wave 
in lead I. Also, the right half of the S matrix, r' and t', 
can be obtained by having an incoming wave in the lead 
II instead of the lead I. But, in our geometries, it can be 
done by changing the direction of the magnetic flux. 

The conductance for this system is obtained from the 
transmission amplitudes, tij, using the Landauer formula. 
For transmission from lead I to lead II the Landauer for- 
mula is given by 
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On the other hand, the conductance from lead II to lead 

I is given by G H ^i = (2e 2 /h) ^5=1 W- But > in tw0 " 
lead conductors, Gz->u should be equal to G//_>j due to 
the reciprocity of S matrix (24). 

IV. Numerical Results 

Let us now show some of our results. The value of A is 
set to 20° throughout, 7 is 180° or 90°. Here e is varied 
from e to e 2 where eo = 0, ei = 40.91, and £2 = 163.7. 
The conductance is calculated using the Landauer for- 
mula, Eq. (22). We define the transmission probability 
Ti = "Y^jZi \tji\ 2 [i < N p ) when there is an zth incoming 
wave in lead I. For each i, Ti + Ri — 1 is satisfied due 
to the unitarity of the S matrix when Ri = \ r ji\ 2 - 
Then the conductance from lead I to lead II is 

2e 2 ^ 

We will compute Ti and T2 for eo < e < £2 for the two 
different lead geometries. 

Figure 3 shows our results for the conductance for the 
case of straight leads for eo < e < £2 when 7 = 180°, 90° 
and a = 0, 0.25. Only Ti is plotted versus e. Note that 
for straight leads, Ti is zero for e < ei, and T2 is zero for 
e < £2- 

Figure 4 shows the analogous results for the case of 
conic leads for e < e < e 2 when 7 = 180°, 90° and 
a = 0, 0.25. We see that Ti (T2) starts to emerge when 
e is lower than t\ (£2)- It is useful to note that this case 
has clear computational advantages. When we integrate 
over Cj's, the Bessel functions are constant, while it is not 
true in the case of straight leads. As a result, the compu- 
tation was much faster. Furthermore, we can keep a^0 
also in leads by choosing suitable bases in leads without 
spending more significant computer time. 

There are several conclusions we can draw from these 
results. First, the shape of the lead matters. For the 
conical leads, the particle can tunnel through the circu- 
lar cavity at energies which forbid any propagation in the 
straight leads. On the other hand, away from tunneling 
energies, both lead geometries give us similar results in 
the most of the range we looked at. 

Second, we can see patterns of fluctuations in all cases 
in Figs. 3 and 4. They have peaks and valleys. (Most of 
valleys go down to zero. We will call them "transmission 
zeros" . We believe that valleys (the transmission zeros, in 
particular) are related to energy eigenvalues of the closed 
billiard. Although there is no exact one-to-one correspon- 
dence, we see loose connections, especially in Figs. 3(b) 
and 4(b). Therefore, these transmission zeros (which are 
reflection peaks) can be treated as resonances with the 
cavity. 

Third, our results depend on the relative position of 
the leads. Although the patterns look similar, we can see 



that the locations of peaks and valleys are different for the 
case 7 = 180° and the case 7 = 90°. Transmission zeros 
for the case 7 = 90° are more likely to match with energy 
eigenvalues of the closed billiard (see Figs. 3(b) and 4(b)). 
This is reasonable because when 7 7^ 180°, electrons are 
more likely to feel the circular billiard. 

Fourth, we can see an indirect result of the AB effect. 
When a = 0.25, we see almost twice as many fluctuations 
as when a = 0. This means that, if we increase a from 
zero, fluctuations will split and be shifted at the same 
time. This is due to the AB effect. 

Let us now look at the conductance from another point 
of view. In Fig. 5, we computed Gi^n of the case of 
conic leads assuming the energy is a complex number. 
We see a series of poles near the real axis. Such poles 
have been seen for simple systems relating to transmis- 
sion zeros Q,|2^],(26). However, we cannot establish an 
exact correspondence because our conductance curves are 
too crowded with fluctuations. One interesting feature, 
however, is the big poles which follow a line across the 
complex plane when 7 = 180° (see Figs. 5(a) and 5(c)). 
When 7 = 90°, Figs. 5(b) and 5(d) show that they become 
less regularly distributed. 

Let us now show our results for the special case of 
a = 0.5. In Fig. 6, we show graphs of Ti and T 2 as a 
function of e for the cases of conic leads when 7 = 180° 
(Fig. 6(a)) and 7 = 90° (Fig. 6(b)). The graphs for the 
cases of straight leads show almost the same results, but 
are not shown here. 

In Fig. 6(b), when 7 = 90°, we see the same type of 
fluctuation pattern that we saw in Figs. 3 and 4. In Fig. 
6(a), however, when 7 = 180°, we see the shut-down of 
the conductance in the lower energy range. This can be 
roughly explained using concepts from the theory of path 
integrals. The transmission amplitude from lead I to lead 
II is approximately proportional to the Green's function 
of Eq. (6) from the point at the center of C\ to the point 
at the center of C2 when A is small. The propagator can 
be obtained using the path integral by summing over all 
paths. These paths can be divided into two groups, upper 
and lower paths. Since this particular case has the reflec- 
tion symmetry along an axis parallel to leads, each lower 
path is the exact mirror image of an upper path, and 
they cancel each other due to the AB effect only when 
a = 0.5. This is destructive quantum interference. If 
the propagator vanishes, the Green's function vanishes, 
too. Therefore, the conductance almost vanishes in this 
case. This is true only when the Fermi wavelength of the 
electron is comparable to the width of the opening in the 
lower energy range. As energy gets bigger, this AB effect 
will be less pronounced because of the finiteness of the 
width of leads. 
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V. Conclusions 

In this idealized model, we have calculated the coonduc- 
tance through the circular AB billiard with two leads us- 
ing the Landauer formula. We have shown that the shape 
and the relative position of leads can play important roles 
in transport properties of the open AB billiard. We have 
observed generic fluctuation patterns, which can be inter- 
preted as resonances. Transmission zeros tend to corre- 
spond well with energy eigenvalues of the closed billiard, 
especially when 7 = 90°. We have also calculated the 
conductance for complex energies, and found regularly 
located big poles. We have seen the significant AB effect 
in the special case of 7 = 180° and a = 0.5. Then the 
conductance was suppressed to zero in the lower part of 
the energy range we looked at. 
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Figure 1: Geometries of two models, (a) Circular AB 
billiard with two straight leads, (b) Circular AB billiard 
with two conic leads. In both cases, the magnetic point 
flux a is placed at the center. R is the radius of the 
billiard, 7 is the angle between two leads, and A is the 
opening angle of leads. 



Figure 2: Energy eigenvalues of the closed circular 
AB billiard as a function of the magnetic flux a(— 
$/<I>o). ein is a dimensionless quantity defined by e;„ = 
(mR 2 /Ti 2 )E ln . 
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Figure 3: T\ (transmission probability for the first mode 
incoming) vs. Fermi energy e for the case of straight leads 
when A = 20°. (a) a = 0, 7 = 180°. (b) a = 0, 7 = 90°. 
(c) a = 0.25, 7 = 180°. (d) a = 0.25, 7 = 90°. Points(o) 
show the energy eigenvalues of the closed AB billiard with 
corresponding a value, e = (mi? 2 /Ti 2 )Ef- e is in the 
range of e < e < £ 2 where e = 0, e\ = 40.91, and 
£2 = 163.7. Ti(i > 1) is equal to zero in this range. The 
conductance Gj^u is obtained from (2e 2 /h)Ti. 



Figure 4: T\ and T 2 (transmission probability for the 
first or second mode incoming, respectively) vs. Fermi 
energy e for the case of conic leads when A = 20°. (a) 
a = 0, 7 = 180°. (b) a = 0, 7 = 90°. (c) a = 0.25, 7 = 
180°. (d) a = 0.25, 7 = 90°. Points(o) show the energy 
eigenvalues of the closed AB billiard with corresponding 
a value. Ti (solid) starts to emerge before ei, and T 2 
(dotted) emerges before £2- The conductance Gi^u is 
obtained from (2e 2 /h)(T 1 + T 2 ). 
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Figure 5: Conductance Gi^n vs. complex energy cr + 
iej for the case of conic leads when A = 20°. (a) a = 0, 
7 = 180°. (b) a = 0, 7 = 90°. (c) a = 0.25, 7 = 180°. 
(d) a — 0.25, 7 = 90°. Dark spots represent poles in the 
complex plane. Small poles near the real axis (ej = 0) 
correspond to fluctuations we saw in Fig. 4. But, we also 
observe big poles marching down regularly. 




Figure 6: T\ and T 2 vs. Fermi energy e for the case of 
conic leads when A = 20° and a = 0.5. (a) 7 = 180°. 
(b) 7 — 90°. Points(o) show the energy eigenvalues of the 
closed AB billiard with a = 0.5. In (a), T x (solid) and T 2 
(dotted) are suppressed to zero in the lower region due to 
the Aharonov-Bohm effect. But, in (b), fluctuations are 
observed as in Fig. 4. 
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